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ABSTRACT 

The binary companion to the peculiar F supergiant HD 172481 is shown to be a Mira 
variable with a pulsation period of 312 days. Its characteristics are within the normal 
range found for solitary Miras of that period, although its pulsation amplitude and 
mass-loss rate Af ~ 3 x 10"^ M^yr"^ are higher than average. Reasons are given for 
suspecting that the F supergiant, which has L ~ 10^ L©, is a white dwarf burning 
hydrogen accreted from its companion. 
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1 INTRODUCTION 

HD 172481 is a high galactic latitude F-type supergiant 
which has drawn some attention since the IRAS survey 
showed it to have an infrared dust excess (Odenwald 1986). 
Reyniers & Van Winckel (2000, henceforth RVW) have 
since shown that the star is a spectroscopic binary with a 
cool M-giant companion. Their detailed spectroscopic anal- 
ysis of the F supergiant indicates a moderate metal defi- 
ciency ([Fe/H]=-0.55), no CNO-enhancement, a slight en- 
hancement of s-process elements and high lithium content 
(loge(Li) = 3.6). Essentially nothing was known about the 
M star prior to the present paper. 

Since the IRAS measurements most authors have as- 
sumed that the F supergiant is a post-AGB star (RVW and 
references therein). This paper first discussed the nature of 
the M giant (section 3) and then, more speculatively, pro- 
vides a possible explanation for the F supergiant (section 4). 



2 PHOTOMETRY 

The published Stromgren and Geneva photometry (Houk 
& Mermilliod 1998) indicate V between 8.9 and 9.1 while 
Geneva photometry from RVW shows that the visual mag- 
nitude is variable with a peak to peak amplitude of 0.22 
mag. RVW show that the M star contributes about 10 per- 
cent of the flux at V when the Mira is at maximum light 
(see below). The value, V — 9.09 mag, given in the SIMBAD 
database is used in the following analysis as representative 
of the mean light of the F star. 

Forty-five JHKL measurements were made between 
1988 and 1999 using the Mkll photometer on the 0.75-m 
telescope at SAAO Sutherland. These are on the SAAO 
system (Carter 1990) and are listed in Table 1. They are 
accurate to about ±0.03 mag at JHK and ±0.05 mag at L. 



3 NATURE OF THE COOL STAR 

A Fourier analysis of the data in Table 1 reveals a clear 
period of 312 days and light curves folded on this period 
are shown in Fig. 1. The peak to peak amplitudes of the 
Fourier fits to these variations are: AJ = 0.53, AH = 0.68, 
AA' = 0.62 and AL = 0.58 mag. These of course are lower 
limits to the amplitudes of the M giant as the F supergiant 
will contribute to the observed magnitude, particularly at 
the shorter wavelengths. Corrected amplitudes are discussed 
in section 3.1 and listed in Table 2. 

These large amplitude variations clearly demonstrate 
that the M giant is a Mira variable and hence at the very 
tip of the Asymptotic Giant Branch (Feast & Whitelock 
1987). 



3.1 Infrared Colours 

The reddening from HD 172481 has been estimated as E{B — 
V) = 0.44 (RVW) and 0.27 (Bersier 1996). In the following 
we calculate the magnitudes and colours using both of these 
values, generally giving the value for E{B — V) = 0.44 first 
and following it with the value associated with E(B — V) = 
0.27, in brackets. 

The first two lines of Table 2 list the magnitudes for 
HD 172481 at the Fourier maximum and minimum (note 
that, as seen in Fig. 1, the phase shifts between the JHKL 
light-curves are negligible). The next three lines have been 
corrected assuming a reddening of E{B — V)= 0.44 and 
give, respectively, the reddened corrected magnitudes of the 
F2Ia star (from Johnson 1966) on the assumption that this 
is the only significant contributor at V, the reddening cor- 
rected mean magnitude of the Mira after the contribution 
from the F star has been removed, and the peak to peak 
amplitude of the Mira variations, also after removing the 
contribution of the F star. These three lines are then re- 
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Table 1. Observed Magnitudes 
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Table 2. Mean Magnitudes 
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Figure 1. Tlie infrared observations for HD 172481 phased on a 
period of 312 days, with an arbitrary zero-point of JD 2440000; 
each point is plotted twice. 



peated using the alternative, E{B ~ V) = 0.27, reddening 
correction. 

It is instructive to compare the colours deduced for 
this Mira with those of normal isolated Miras. Whitelock, 
Marang & Feast (2000, henceforth WMF) provide JHKL 
observations of Miras with thin dust shells, while Whitelock 
et al. (1994) do the same for IRAS selected Miras with rel- 
atively thick shells. The mean colours of the Mira are listed 
in Table 3 together with those of a typical 312 day Mira, 
calculated from equations 1 to 4 of WMF. 

The deduced colours are obviously somewhat uncertain, 
particularly at J, because of the uncertain correction for the 
F star. Nevertheless, they fall within the range shown by nor- 
mal unreddened Miras. In particular, K — L, although larger 
than the mean calculated from WMF equation 3, is within 
the range shown by normal Miras. WMF note a strong cor- 
relation between the pulsation amplitude and K — L, and 
indeed the pulsation amplitude of this Mira is high. 

This Mira is of particular interest in view of the fact 
that we can infer its metallicity; only a very small number 
of Miras, those associated with globular clusters, have reli- 
able metallicities. RVW estimated [Fe/H]=-0.55 for the F 
supergiant which we might reasonably assume also applies to 
the Mira. For comparison. Table 3 gives the colours (from 
WMF) of the Mira V3 in NGC5927 which has a 311 day 
period and [Fe/H]=-0.37 (Harris 1996). 
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Table 3. Mean Colours of various components 
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3.2 Mass-loss Rate 

As RVW discussed, HD 172481 was detected by IRAS and 
shows an obvious infrared excess. The colour corrected 
Ufim mag is [12] = 1.76, so K - [12] = 4.3 (4.2). This is 
larger than values shown by any of the stars discussed by 
WMF, but is within the range found for IRAS Miras by 
Whitelock et al. (1994). This K - [12] would imply a mass- 
loss rate of M ~ 3 x 10"® Moyr"^ (Whitelock et al. 1994 fig 
21) which is high for a 312 day Mira, but higher mass-loss 
rates are usually associated with larger pulsation amplitudes 
and again we note that the pulsation amplitude of this star 
is high. 

It would therefore seem possible that the dust shell, 
which is responsible for large A!"— [12] and presumably for the 
circumstellar extinction, originates entirely from the Mira. 
While it is not essential to assume that any of the dust is 
associated with the F supergiant, we cannot rule out the 
possibility that some of it is. Neither can we rule out the 
possibility of dust trapped in a circumbinary disk, as has 
been proposed for some other supposed post-AGB binary 
systems (e.g. Waters et al. 1993; Van Winckel 1999) 

Note also that the IRAS colours of HD 172481 are 
within the range found for D-type symbiotics (Whitelock 
1988), illustrating the similarity of their dust shells (see sec- 
tion 4). 

3.3 Distance and Luminosity 

The distance modulus of the Mira can be calculated via the 
period luminosity relation (Whitelock & Feast 2000): 

Mk = -3.47 log P + 0.84, 

to give {Ko - Mk) = 13.9 (13.8) or a distance of 6.0 (5.7) 
kpc, and a height above the galactic plane of -1.1 (-1.0) kpc. 

The bolometric magnitude of the Mira can be es- 
tablished by fitting a blackbody to the JHKL flux (see 
Robertson & Feast 1981). This gives a bolometric magni- 
tude at mean light of 9.30 (9.09) and at maximum light 
of 8.72 (8.64). Or in absolute terms at mean light Mboi = 
-4.59 (-4.69) and Lm = 5.3 x 10^ (5.8 x 1O^)L0. 



4 THE F SUPERGIANT 

RVW established that the relative luminosities of the F su- 
pergiant and the M star are Lf/Lm = 1.8, for a reddening 
of E{B ~V)= 0.44 (or Lf/Lm = 1.0 for E{B -V)= 0.2). 
The near-infrared measurement which they used to de- 
fine the fiux of the system was obtained on JD 2448850 
[K = 5.48 Van Winckel, private communication), i.e. at 
maximum light for the Mira. Thus the luminosity of the F 
star is Lf = 1.6 x 10'' Lq (0.9 x 10'' Lq). 

RVW argue that the F supergiant is a post-AGB star 



on the basis of its high galactic latitude, large radial velocity 
and moderate metal deficiency. In support of this they also 
mention the dust, a variable Ha profile and a slight over- 
abundance of s-process elements. Our identification of the 
M star as a Mira variable confirms the binary as a low mass 
system. There is, however, an alternative to the post-AGB 
star scenario. The F supergiant could be a white dwarf, with 
hydrogen-shell burning of material accreted from the Mira 
wind providing a pseudo-photosphere and the observed lu- 
minosity. 

Iben & Tutukov (1996) discuss the evolution of mass ac- 
creting white dwarfs in wide binaries systems, such as sym- 
biotic stars. Of particular interest in the present context are 
the D-type symbiotics where a white dwarf accretes from the 
wind of a Mira variable in a binary system with a period in 
excess of about 10 years (Whitelock 1987). In certain cases 
the accretion rates are such that a thermonuclear runaway 
gives rise to a very slow nova, e.g. RR Tel and V1016 Cyg. 
These symbiotic stars, with their high excitation emission 
lines, are apparently very different from HD 172481. How- 
ever, it is clear from Iben & Tutukov that relatively small 
differences, e.g. in the accretion rate of the white dwarf, can 
give rise to very different phenomena. 

The luminosity derived above, Lf ~ 10* Lq, is that 
expected for a cool ~ 0.6 Mq white dwarf in a hydrogen shell 
burning phase following a hydrogen shell flash (see Iben & 
Tutukov fig. 5). 

A very rough estimate of the accretion rate of the white 
dwarf can be derived from Iben & Tutukov's equation 2. If 
we assume the mass-loss rate from the Mira (section 3.2), a 
0.6 Mq white dwarf, a total stellar mass (white dwarf plus 
Mira) of 1.5 Mq, a wind velocity of 20 kms"'^ and an orbital 
period of 10 years, the derived accretion rate onto the white 
dwarf is 2 X lO"^M0yr"\ 

It is clear from fig. 7 of Iben & Tutukov that this rate 
would in fact result in steady hydrogen burning on the sur- 
face of the white dwarf. However, for the parameters as- 
sumed this rate should probably be regarded as an upper 
limit, as Iben & Tutukov's equation 2 is based on Bondi- 
Hoyle accretion and other estimates lead to significantly (up 
to ten times) lower values (e.g. Theuns, Boffin & Jorissen 
1996; Mastrodemos & Morris 1999). Furthermore, plausible 
changes in the period or the wind velocity could increase or 
decrease the rate. Nevertheless it seems likely that a white 
dwarf paired with this Mira would accrete mass in such a 
way as to either burn it steadily or in shell flashes, for a 
large range of stellar separations. 

In discussing the nova phenomenon Iben & Tutukov 
point out that in wide binaries there is no reason to expect 
the loss of the hydrogen envelope after the explosion. So 
that it is possible for the newly created supergiant to last 
for several decades. This might explain the present condition 
of the F supergiant in HD 172481. 

The high abundance of lithium in the supergiant spec- 
trum (RVW) may well be key to understanding this sys- 
tem. Arnould & N0rgaard (1975) and Starrfield et al. (1978) 
predict the formation of large quantities of lithium during 
hydrogen-burning themonuclear runaways. However, their 
calculations are not obviously applicable to the situation 
under discussion and they also predict an overabundance 
of carbon and nitrogen which are not found in HD 172481 
(RVW). More recently Hernanz et al. (1996) predict signif- 
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icant lithium production from some nova models, but these 
seem to underestimate the ejecta masses. Clearly more the- 
oretical work is required in this area. 



ACKNOWLEDGMENTS 

Our thanks to Hans Van Winckel for providing a preprint 
of their paper and details of their unpublished infrared pho- 
tometry. We are also grateful to Michael Feast for helpful 
discussions. This paper is based on observations made from 
the South African Astronomical Observatory and also made 
use of information from the SIMBAD data-base. The com- 
ments and suggestions of an anonymous referee were also 
helpful. 



REFERENCES 

Arnould M., N0rgaard H., 1975, AA, 42,55 
Carter B. S. C, 1990, MNRAS, 242, 1 
Bersier D., 1996, A&A, 308, 514 

Feast M. W., Whitelock P. A., 1987, in: Late Stages of Stellar 
Evolution, (eds.) S. Kwok, S. R. Pottasch, Reidel, Dordrecht, 
p. 33 

Har ris W. E.. 1996. AJ. 112. 1487 

(http: / /physun. physics. mcmaster.ca/Globular.htm^ ) 
Hauk B., Mermilliod M., 1998, A&AS, 129, 431 
Hernanz M., Jose J., Coc A., Isern J., 1996, ApJ, 465, L27 
Iben I., Tutukov A. V., 1996, ApJS, 105, 145 
Johnson H. L., 1966, ARAA, 4, 193 
Mastrodemos N., Morris M., 1999, ApJ, 523, 357 
Odenwald S. F., 1986, ApJ, 307, 711 
Re yniers M.. Van Wi nckel H., 2000, A&A, in press, 

astro- ph/0010486| (RVW) 
Robertson B. S. C, Feast M. W., 1981, MNRAS, 196, 111 
Starrfield S., Truran J. W., Sparks W. M., Arnould M., 1978, 
ApJ, 222, 600 

Theuns T., Boffin H. M. J., Jorissen A., 1996, MNRAS, 280, 1264 
Van Winckel H., 1999, in: Asymptotic Giant Branch Star, lAU 

Symp. 191, p. 465 
Waters L. B. F. M., Waelkens C, Mayor M., Trams N. R., 1993, 

A&A, 269, 242 
Whitelock P. A., 1987, PASP, 99, 573 

Whitelock P. A., 1988, in: The Symbiotic Phenomenon, (eds.) J. 

Mikolajewska et al., Kluwer, Dordrecht, p. 47 
Whitelock P. A., Feast M. W., 2000, MNRAS, 319, 759 
Whitelock P. A., Menzies J. W., Feast M. W., Marang F., Carter 

B., Roberts G., Catchpole R. M., Chapman J., 1994, MNRAS, 

267, 711 

Whitelock P. A., Marang F., Feast M. W., 2000, MNRAS, 319, 
728 (WMF) 



© 2001 RAS, MNRAS 000, 



